Quark number susceptibilities, strangeness and dynamical confinement 
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We report first results on the strange quark number susceptibility, Xs , over a large range of tem- 
peratures, mainly in the plasma phase of QCD. Xs jumps across the phase transition temperature, 
Tc, and grows rapidly with temperature above but close to Tc. For all quark masses and suscep- 
tibilities in the entire temperature range studied, we found significant departures from ideal-gas 
values. We also observed a strong correlation between these quantities and the susceptibility in the 
scalar/pseudo-scalar channel, supporting ideas of "dynamical confinement" in the high temperature 
phase of the QCD plasma. 
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As seen recently in the Quark Matter 2001 conference 
1^ , the first runs of the RHIC have provided many excit- 
ing results. With the anticipated much longer runtime, 
at a variety of colliding energies and nuclear species, the 
existence and properties of the high temperature phase 
of quantum chromodynamics (QCD) will be probed and 
tested further. This is an immediate reason to improve 
our theoretical knowledge of such matter. 

Fully non-perturbative computations of quark number 
susceptibilities are important for three reasons. Firstly, 
there have been attempts to link them directly to exper- 
imental measurements of event-to-event fluctuations in 
particle production Q. Secondly, experimental observa- 
tions of a relative enhancement of strange quarks have 
been attributed to the formation of a QCD plasma ||]; a 
hypothesis which can be quantitatively tested against the 
computation of the strange quark susceptibility. Finally, 
earlier results |^-^ showed a strong jump in the suscep- 
tibility across the phase transition, but indicated a sta- 
tistically significant 20% departure from weak-coupling 
behaviour; this is the only (lattice) evidence that physics 
at finite chemical potential is not weak-coupling physics. 
In this paper we present the first extensive systematic 
study of such susceptibilities and relate them to other 
widely studied variables. 

The partition function of QCD with three quark 
species can be written on a discrete space-time lattice 
as 



exp[-S'(r)] X 
n det M(r,m/,/i/), (1) 



where the temperature T enters through the size of the 
Euclidean time direction, S{T) is the gluonic part of the 
action, the determinants of the Dirac matrices, M , con- 



tain as parameters the quark masses, m, and the chem- 
ical potentials /i. The chemical potentials for specific 
flavours, /i/ (/ = m, d, s), can be expressed as linear com- 
binations of the chemical potentials, (a=0,3,8), corre- 
sponding to the diagonal generators of flavour SU{3) . We 
work with the choice Ho = iJ-u + l^'d + IJ-s, Ii-s = fJ-u — l^d 
and /ig = + fJ-d — '^Ii's, because hq is then the usual 
baryon chemical potential. 

We follow convention in defining quark numbers 



and number susceptibilities 
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Here V denotes the spatial volume, and i and j are ei- 
ther of the index sets / or a. Transformations between 
these bases are straightforward. We lighten our notation 
by writing the diagonal susceptibilities xu as Xi- We de- 
termine the susceptibilities at zero chemical potentials, 
fif = 0. in this limit, of course, ni{T) = for all i. It 
is also easy to prove that xo3 = if rn„ = rrid . Since all 
current lattice computations are made in this approxima- 
tion, the off-diagonal susceptibility cannot be measured 
at present. Also, in this limit, Xus = Xds, where the 
strange quark mass, rus ^ rriu^d- 

Previous works have reported the values of the flavour 
SU{2) singlet, xoj and triplet, X3, susceptibilities. Mea- 
surements of xo,3 in 2-flavour dynamical QCD showed a 
statistically significant 20% departure from free field the- 
ory (FFT) for m/T = 0.1 and 0.15 Later compu- 
tations with m/T = 0.05 reported xo 3 compatible with 
FFT but with errors of about 20% Quenched com- 
putations deviated from the first set of dynamical QCD 
results by 5-10% [Q. All this work covered a range of 
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temperatures up to l.bTc- However, the quark mass var- 
ied with temperature, since the quantity that was fixed 
was m/T. 

In order to faciUtate a comparison of our results with 
these, we report the same susceptibihties. In the re- 
mainder of this paper, xo and X3 refer to these 2-flavour 
quantities. We also report the first determination of the 
strange quark susceptibility, Xs- This work improves on 
previous studies in four other ways — by covering a larger 
range of temperatures, by using a series of different quark 
masses at fixed m/Tc, by investigating finite spatial vol- 
ume effects systematically, and by analysing the effect 
of increasing statistics in the stochastic determination of 
Fermion operators on the lattice. The study of volume 
dependence gives control over extrapolation to the ther- 
modynamic limit. Our study of many different quark 
masses gives the strange quark susceptibility. The sys- 
tematic study of the stochastic method yields a vastly 
improved determination of the flavour-singlet suscepti- 
bility, Xo- 

These computations have been made on lattices with 
lattice spacing a — 1/ 4T in the quenched approximation. 
Fermion loops are neglected in this approximation, mak- 
ing it substantially easier to handle numerically. Apart 
from an overall normalisation of the temperature scale, 
this approximation is known to reproduce all the quali- 
tative features of the full QCD simulations. Furthermore 
numerical agreement between the quenched and 2-flavour 
dynamical QCD results for xo,3 are obtained by 5-10% 
correction of the former It has been shown recently 
that one can extract continuum results from the lattice 
spacing we employ Nevertheless, in future studies 

we will analyse the effects of relaxing these two approxi- 
mations. 

Successive configurations used in our computations are 
separated by 1000 sweeps of a Cabbibo-Marinari heat 
bath algorithm, so that the gauge fields are completely 
decorrelated by any measure one may choose to use. At 
P = 5.8941, corresponding to T = 1.5Tc we have gener- 
ated configurations on 4 x S'^, 4 x 12^^ and 4 x 16'^ lat- 
tices. At other couplings, corresponding to T = 0.75Tc, 
l.lTc, 1.25Tc, 2Tc and 3Tc, we have used only the 4 x 12^ 
lattice. We have used quark masses m/Tc = 3, 1.5, 1, 
0.75, 0.30 and 0.03. Using estimates of re=275-290 for 
the quenched theory [|o|, we see that the strange quark 
mass lies between 0.3Tc and 0.5rc. 

In this letter we use staggered Fermions exclusively. 
Since these are defined for four flavours, we have to nor- 
malise the susceptibilities appropriately pl|. We have 
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where the two operators involved are 



T 



d = — (Tr {M"Ivr^ - M'M-^M'M-^)) 
O2 = ^((TrM'M-i)' 



(5) 



The traces here are sums over lattice points and colour in- 
dices, and angular brackets are averages over gauge field 
configurations. Primes denote derivatives of the Dirac 
matrix with respect to appropriate chemical potentials. 
The quark mass to be used in the Dirac operator for 
evaluating Xs is, of course, different from that for Xo,3- 

The traces are evaluated by the usual stochastic tech- 
nique. 
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where Ri are a set of N uncorrelated vectors with com- 
ponents drawn independently from a Gaussian ensemble 
with unit variance. Each vector has three colour com- 
ponents at each site of the lattice. We improve on the 
definitions in eq. (||) by using half lattice versions of the 
Dirac operator for staggered Fermions. A detailed dis- 
cussion of the stochastic evaluation of the squared trace 
in eq. (||) can be found in . A systematic study of the 
optimum value of N is shown in Figure |l]. We find that 
« 80 is needed in order to see that xo = Xs within 
statistical errors jl^] . In all our subsequent work we have 
used A^ = 80 ^ . Such a large value of A^ also seems to 
decrease the variance in the average over gauge configu- 
rations. 
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FIG. 1. Variation of xo (circles) and X3 (squares) with the 
number of vectors, A'^, used per configuration. The averages 
are taken over 20 thermal configurations at T = l.STc on a 
4x8^ lattice for m/T^ = 1. 

In Figure |^ we exhibit the spatial volume dependence 
of our results at T = 1.5Tc. Note that the volume de- 
pendence is smaller than the statistical errors on x at 
this temperature, and also much smaller than what could 
be expected for an ideal gas of Fermions. In view of 
this, we have chosen to perform the remaining compu- 
tations with lattices of size 4 x 12'^. This volume allows 



2 



us to measure thermodynamic quantities up to 3Tc. At 
larger T , spatial deconfinement sets in and distorts the 
results unless larger lattices are used . Measurements 
of finite volume effects on the QCD equation of state in- 
dicates that this lattice size can also be used down to 
l.lTc, below which the first order phase transition of the 
quenched theory causes strong finite volume effects. As a 
result, the quenched theory approximation to full QCD 
is expected to fail close to T^- Another finite size ef- 
fect appears in Fermion computations in the quenched 
theory. As the quark mass decreases at fixed tempera- 
ture, the scalar/pseudo-scalar screening length increases. 
If this length becomes comparable to the spatial dimen- 
sions of the lattice, then finite volume effects cannot be 
neglected. Among all our computations, this happened 
only for m/T^ = 0.03 at T l.lTc- 
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FIG. 2. Quark number susceptibilities at T = 1.5Tc. Val- 
ues of xo are displaced slightly to the right for visibility. 

In Figure |^ we collect together all our data. At all 
these points, Xo = X3 = "^Xs within statistical errors Q . 
For T < Tc all the susceptibilities are consistent with 
zero. Notice that the results lie significantly below the 
expectation from FFT even at temperatures as high as 
3Tc. It is interesting to note that the departures from an 
ideal gas become stronger with increasing quark mass. 
We discuss later our checks that this is not a lattice ar- 
tifact. The jump in Xs across Tc, its rapid increase with 
temperature, and its becoming comparable with xo,3 for 
T > 2Tc, all have observable consequences in the pattern 
of relative strangeness enhancement, including strange 
baryon enhancement, in going from CERN to RHIC en- 
ergies ll5| . 

It is interesting to speculate upon the reasons for 
departure from the weak coupling theory. In view of 
the long-standing observation that screening masses in 
the scalar/pseudo-scalar channel (so-called pion screen- 
ing masses) also depart strongly from the weak coupling 
theory, we would like to check whether these observations 
are related. Below Tc, at vanishing chemical potential, 
pair production of quarks can take place only by pair 
production of mesons. The lightest meson, the pion, will 
be most effective at producing quark pairs. Hence the 
pion susceptibility 



G,(fco-O,k = 0) = ^^a(z) (7) 



(here is the pion propagator in momentum space and 
is the zero-momentum screening correlator) should 
determine xo,3- For T > Tc, this logic does not neces- 
sarily follow, but in view of the strong correlation in the 
pion sector, it is interesting to test such a hypothesis. 




FIG. 3. The ratio of xa to its value in a free fermion theory 
at the same lattice volume and quark mass, shown as a func- 
tion of temperature for the values of mq/Tc indicated. The 
lines are cubic spline fits to the data. The shaded area is the 
range of masses relevant for Xs- 

All the susceptibilities are functions of m and T. Since 
QCD with two light dynamical flavours is expected to 
have a second order chiral phase transition for m — Q 
and T — Tc, near the critical point we should expect these 
functions to scale as powers of T — Tc or m. These pow- 
ers are determined by the ratio of various critical expo- 
nents. The scaled quantities would then vary as universal 
functions of each other. At the critical point, there is a 
zero mass correlation in the scalar / pseudoscalar channel 
which is solely responsible for the susceptibilities, since 
the correlations in the other channels are massive. The 
universal scaling is then a signal for a restricted form of 
dynamical confinement near the critical point. 

In quenched QCD, where there is no critical point, we 
can instead make an expansion such as 

Xi = a°{T) + a\{T)m + aj{T)m\ogm + 0{m'^), (8) 

where i = 0, 3, s or tt |1^ ]. Eliminating m between Xn- 
and Xs, we can always expand one in terms of the other, 
albeit with temperature dependent coefficients. The sur- 
prise, shown in Figure ^, is that the curves for different 
T can be scaled to lie on top of each other. The scal- 
ing function which achieves this goes to 1 as T —> Tc. 
Thus, for T < 2Tc and m/T > 1.5, the co-variance of 
Xtt and X3 at Tc determines that at higher temperatures. 
Since the scalar/pseudoscalar screening mass is smallest 
at finite temperature, and very much smaller than other 
screening masses, therefore the observed scaling may be 
interpreted as a demonstration of dynamical confinement 
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in much the same way as at a second order phase transi- 
tion. 

Why is this not a trivial observation? After all, as 
~* 0, Xir/X3 is a T-dependent number, which defines 
the scaling function. The point is that the scaling is 
observed not only for m = 0, but for a whole range of 
masses. This defines an universal curve at Tc, not just 
a single value. As a result, this scaling also gives us a 
partial ability to test for lattice artifacts. The very fact 
that scaling is seen for all m/Tc > 1.5, when results for 
the lightest mass are consistent with previous computa- 
tions in two flavour dynamical QCD, implies that there 
are neither strong quenching artifacts, nor other strong 
lattice artifacts in our measurements. 
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FIG. 4. Scaled values of Xtt vary with xs independently of 
the values of T and m/T at which they are measured, pro- 
vided T < 2Tc and m/Tc < 1.5. The meanings of the symbols 
are the same in the two panels. X3 increases with decreasing 
m/Tc at fixed T. 

In summary, we have presented new and precise re- 
sults on quark number susceptibilities over a wide range 
of temperatures and quark masses in the high temper- 
ature phase of QCD. The susceptibilities differ signifi- 
cantly from the ideal gas expectations (Figure ||) . These 
deviations increase with mass and decrease at higher T. 
As a result, we expect the relative strangeness enhance- 
ment seen in heavy-ion collisions to increase with tem- 
perature as quantitatively determined here. The linear 
relation between x-k a-nd xs (Figure ^) is the clearest ev- 
idence to date for the hypothesis of "dynamical confine- 
ment" in the high temperature phase of the plasma . 
However, it also shows that such a picture becomes less 
effective with increasing temperature. It is perhaps not a 
coincidence that the temperature at which the quenched 
plasma becomes free of this phenomenon is also the tem- 
perature at which dimensional reduction becomes quanti- 
tatively correct For T > 2Tc, the light quark suscep- 
tibility is known accurately enough to test dimensional 
reduction and various other ideas which have emerged in 
trying to explain lattice results on the QCD equation of 
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